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TECHNICAL PAPER 


INHERENT CONSERVATISM IN DETERMINISTIC 
QUASI-STATIC STRUCTURAL ANALYSES 


I. INTRODUCTION 


In designing reliable and affordable next generation carriers to orbit, performance of aerospace 
frame structures once more emerges as a critical limiting component requiring innovative configurations, 
updated materials, and refined analyses. And because designing for structural safety is in direct contention 
with performance and cost, a study was initiated to explore the prevailing deterministic structural safety 
factor assumptions and standards for their often suspected conservatism. 

The origin of the conventional safety factor, its simplistic application, its verification criteria, and 
launch vehicle performance sensitivities to its excesses were briefly reflected upon. Though raw data 
dispersions were noted to be processed through probability techniques, current practice is to reduce them 
to deterministic input values for quasi-static substructural loads and stress computations; and there lies the 
cradle of excessive conservatism. As deterministic values were combined and manipulated throughout the 
computations, the imbedded statistical dispersions were impelled to be summed rather than root-summed- 
squared, which violated error propagation laws that predicted cumulative excessively applied loads and 
stresses. 

Conserving the data in statistical format and complying with the error propagation laws throughout 
the structural computation processes lead to designer controlled and leaner safety factors which are treated 
in the conventional manner to assess the structural relative safety and to experimentally verify the structural 
response. In examining the role of the safety factor in the failure concept, the safety factor only increases 
the number of standard deviations of the applied stress, which decreases the applied and resistive stress 
overlap. 

The safety factor expressing the applied and resistive stresses in statistical format was integrated 
into the first order reliability method to provide the option of designing the structural frame to a specified 
uniform and absolute reliability, or to relate the arbitrarily specified design safety factor to a normalize 
reliability. Unexpectedly, the safety factor as currently applied is proportionately relative only with materials 
and welds having identical coefficients of variation. As the coefficient of variation increases, the relative 
proportion of the safety factor and the normalized reliability decreases. 



II. CONVENTIONAL SAFETY FACTOR 


The conventional safety factor has served the aerostructural community very well through many 
progressive changes in materials, associated disciplines, and subjective assumptions. Perhaps its success is 
owed to its unchallenged safety conservatism and to its neglected performance potential. It should be 
instructive to assess excessive safety factor consequences to payload performance, and subsequently to 
cost, through its sensitivity to surplus weight over global structural areas. Local stress concentration regions 
identified by abrupt changes of geometry, loads, metallurgy, and temperature are of less consequence to 
performance. 


A. Historical Note 

Though there was very little statistical data on loading conditions in 1932, there was evidence that 
successfully designed airplanes did not yield. Since the common structural materials were 1 7ST aluminum 
alloy and 4130 steel, and since the aluminum had an ultimate-to-yield stress ratio of 1 .5, the arbitrary 1 .5 
safety factor at fracture was universally accepted 1 by the Commerce, Air Force, and Navy departments. 
Steel structures were obviously penalized by that standard, but it was acceptable because of their limited 
applications at that time. During the Apollo program, NASA field centers chose to reduce the ultimate 
safety factor to 1.4, in order to capitalize on recently improved aluminum properties to increase vehicle 
performance. Again its penalty to high strength steels was ignored. The 1 .4 factor at fracture is now an 
official NASA standard, 2 and it is expressed by: 


SF =1 A-— 
* r tu s A 


( 1 ) 


Though safety factors generally are specified at all levels of material fundamental property changes, 
the safety factor based on polycrystalline yield is difficult to verify. Plastic deformation starts in different 
locations, numbers, and intensities, and it is hard to detect and determine where and how much deformation 
has progressed until large enough parts have been affected and detected. This phenomenon explains why 
different gauge lengths in uniaxial tensile tests provide different elastic limits, why yield coefficients of 
variations are higher than strength variations, and why the elastic limit is more difficult to detect in brittle 
materials. Exceeding the yield point permanently changes the structural boundary conditions and reduces 
fatigue life. Similar degradation phenomena may exist in composites (matrix cracking) with varying 
consequence in stiffness or utility. These, among others, are reasons for which explicit yield safety factors 
should be contingent on the consequence-^ of each operational case. 

Nevertheless, current deterministic experimental tests consist of a static structural elastic response 
verification to predicted maximum operational environments, and of an ultimate safety factor of 1.4 
verification to avoid operating in the plastic region of most aerostructural polycrystalline materials. The 
fracture safety factor subsequently has been rationalized to cover rare operational events in which no 
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statistical design data exists. Its traditional and historical usage now exerts the greatest influence on design 
and contractually binding acceptance criteria. 

It should be cautioned that a deterministic static test is a pass-fail experiment of combined physical 
phenomena resulting in single numbers that often approximate expected fracture values for the wrong 
reasons. Comprehensive pretest and posttest analyses are essential for resolving the right reason and for 
learning the most from the test investment. A test not conducted to fracture provides little information no 
matter how well the complex structure performs thereafter. While the loading instant, location, and nature 
of yielding may be difficult to experimentally detect, testing to fracture leaves little doubt. 

B. Consequence of Excessive Conservatism 

A deterministic static test should prove the article to not be marginally or excessively safe for all the 
right reasons. A meaningful stress audit should present negative and positive margins and consequences 
for both cases. While the cause and consequences of negative margins derived from tests are invariably 
modified, sources and consequences of excessive positive margins are often ignored at the ultimate expense 
of payload performance. The dilemma and magnitude of excessive conservatism may be best appreciated 
through illustrations of the carrier performance sensitivities to surplus structural safety. For example, a 
monocoque shell structural weight is defined by: 

W s = C x rlt 

where two of the dimensions (r and /) envelop the structural element size and shape that are usually optimized 
by the system’s operational environments, payload performance, and cost. The thickness, t, is controlled 
by normal stress limitations and by the designer’s selected safety. Applying the safety factor on the load, p, 
a hypothetical tank minimum thickness is approximated from the material physics: 

f= p(SF)r 

S tu ’ 

and the shell weight is rewritten to relate to the safety factor by: 

C, lr 2 p(SF) 

W s~ S 
°tu 

The weight performance sensitivity 4 to the safety factor is given by the change in weight to change in the 
safety factor: 

3W S _ C,lr 2 pd(SF)S tu _ ;j(SF) 

C, lr*p(SF)S, u <SF) ' 

resulting in a direct proportionality of 1 percent increase in weight for each percent increase in safety 
factor. This sensitivity may be a useful thumb rule for assessing the safety factor penalty to structural 
element performance subjected to normal stresses. 



The length and width of a plate in bending is another example of the element size that is fixed by 
the structural system, and the thickness parameter is approximated from: 


t= 


C 2 M(SF[ 


bS, 


til 


where the ultimate tension and compression stresses are assumed symmetrical. The plate weight is: 


W=CM 


M(SF) 

bS tu 


1 

2 


Proceeding as before, the bending plate weight sensitivity to the safety factor is the partial of the weight 
divided by the weight: 


w s _u j (SF) ’ 

and the bending sensitivity turns out to be half of the normal stress sensitivity. However, bending stresses 
are primarily local and are not as weight prevalent as normal stresses in aerospace structures. Sensitivities 
should be verified on all global, high-performance flight structures to effectively support trades and design 
insight. 


The ultimate ripple effect of excessive safety factors may be realized from flight performance 
parameters. Using the well-known rocket equation: 


W 

A V- A Vloss = I sp gln jy w > 

s p PL 

and assuming the orbital and propulsion parameters are constant, then the mass fraction remains a constant: 


AV-AVloss 
exp. V? 


VV 


— Ca , 


W P +W S+ W PL ^ 


and the propellant weight to orbit is: 


W P = -^< W s + W PL> 
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The sensitivity of the propellant weight increase to accommodate structural weight increase is: 

dW p _ w s dw s 
w p ~w s+ w PL w s ■ 

Using the weight-to-safety factor relationship developed above, the sensitivity of increased propellant 
weight consumption to accommodate the safety factor increase is: 

dW P _ W s d(SF) 

% W S + W PL (S f ) ' 

The ripple effect continues in that, increasing the propellant weight further increases the tank size 
and tank weight, which necessitates more propellant weight, etc. The increased tank size and associated 
propellant loading facilities represent the initial manufacturing costs. The increased tank and propellant 
weights are the recurring costs of lost performance. Recognizing the penalties of excessive safety factors 
and potential rippling effects, then it seems not enough for a senior structural analyst to design a reliable 
structure. His hallmark should be a lean reliable design such as to create and shift the least excessive 
conservatism burden downstream onto the vehicle performance and supporting disciplines. 

On the other side of conservatism is the acceptance and compensation of marginal assumptions 
made in structural modeling. The above shell thickness example was a strength of materials approximation 
of an elasticity theory tube in which the thickness-to-radius ratio is assumed very small, such that the radial 
strain may be assumed uniform across the thickness to simplify the compatibility and boundary conditions. 
The error is often less than 1 percent for common shell properties and is usually ignored. The plate is also 
a strength of materials approximation which assumes the elastic cross section to remain rectangular after 
bending. Most analytical and finite element method models are also approximations to simplify and expedite 
solutions with negligible resulting errors. Nevertheless, it is incumbent on the analyst to assess and estimate 
errors before dismissing them or lavishing approximate solutions with arbitrary gut compensating factors. 
It’s all part of the professional bag of creditable experience. 
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m. DATA CHARACTERIZATION 


Models are idealized into the simplest mathematical expressions within the physical phenomena of 
the data and its intended application. Not all data is equally important, as noted by Pareto’s principal 5 and 
as can be distinguished by sensitivity analyses. Data having negligible effect on performance may be 
reduced to deterministic values. Data having major consequences must be characterized and processed in 
statistical format throughout the computational process involving single and combined phenomena. 

Data that must be expressed and processed in statistical format is developed in this section, and 
those familiar with probability methods may wish to defer it. Those pragmatists interested in reviewing 
basic probability expressions will find most of the statistics and probabilistic techniques required for 
structures in this section . 5 While acknowledging the immense contributions of statisticians to structural 
analyses, it is often more important, in an established process, that structural analysts learn a little statistics 
than a statistician know a little structures. 

A. Central Moments and Distributions 

Variations are intrinsic in all observed phenomena and are of little engineering information in raw 
form. The best approach to summarizing a table of raw data of any distribution is to define the centroid 
about which the data is scattered. This variable is the first central moment 6 of the independent variables 
commonly known as the sample mean, or sample average, and is defined by: 



where x, is the ith specimen value, and n is the total number of specimens. The sample mean is calculated 
from a limited sample size and is, therefore, an estimate of the population mean. A measure of the dispersion 
of the data about the mean is the second central moment known as the sample variance, and its square root: 


<y-\ 


r J&rv) 

1 i=l 


( 3 ) 


is called the sample standard deviation “cr“ and is a measure of the actual variation in a set of data. 

The coefficient of variation is a relative variation, or scatter, among sets and is defined as the ratio 
of the standard deviation and the mean: 



( 4 ) 
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The coefficient of variation is an effective technique for supporting judgment through comparison with 
other known events. Coefficients of variation are known to be small for biological phenomena, but they are 
large for natural materials. Coefficients of variation are small for highly controlled manmade materials and 
are larger for brittle materials. A knowledge of typical coefficients of recurring sources may serve as a 
source forjudging quality and acceptability of data. Estimates of some common material structural properties 
characterized by coefficients of variations are metal ultimate strengths 0.05, yield 0.07, weights 0.015, 
steel fracture toughness 0.07, constant amplitude fatigue 0.40, and filament composite strength 0. 12. 

Another technique used to evaluate raw data is the population probability density distribution, 
which defines the area shape of the distribution. Distribution shapes are fixed by their natural scatter of 
data about their means. Shapes are modeled by distribution functions to estimate the probability of a desired 
value for an assigned range of probability. As shapes become more complex, distribution models become 
more difficult, and skills and labor to apply them escalate. There are obviously many continuous distributions 
that may be constructed and favored by academicians, but normal distributions are most widely used by 
practitioners because they are the easiest and because the mean of independent observations is believed 
to approach a normal distribution as “n” approaches infinity (central limit theory). It is also a good 
representation of many natural physical variables or for small samples with no dominating variance. The 
equation of the normal probability density distribution is: 


f(x)= 




2 

, for — °°<jt<oo , 


(5) 


assuming true values of a very large sample size. Normally distributed phenomena are sometimes disguised 
as nonnormal when data samples are selected from casually broadened and unscreened sources. Most 
metallic mechanical properties are known to be normally distributed, though fatigue properties are not, as 
fatigue is presently understood. 

Analytical advantages in using normal probability distributions are that they may be completely 
defined by two variables (fl and <f), and it is the most developed theory having many of their characteristics 
well established. The area within a specified number of standard deviations of a probability density plot 
represents the proportion of the data population captured. One standard deviation about the average of a 
normal distribution is calculated to capture 68.3 percent of the data tabulated as illustrated in figure 1 . Two 
standard deviations include 95.5 percent of the data, and three standards include 99.7 percent. 



FIGURE 1 . — Normal probability density plot. 
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The normal distribution plot and the mathematical test for the normality of data distribution are 
programmed in quick basic in appendices A 1 and A2 for convenience. Most engineering data distributions 
are one-sided, occurring in the lower or upper sides of the data, such that either side may be developed into 
a split normal, should the distribution not pass the normality test. This cognizance and the central limit 
theory justify the presumption that structural probability demand and capability distributions may be 
normalized by constructing a mirror image of the engaged side about its peak frequency value and calculating 
the standard deviation from the constructed symmetrical distribution. This universal normalizing technique 
is illustrated in appendix A3 and is applicable to normal and nonnormal observed structural data. Hence, all 
structural data in this study is assumed to be generalized into normal probability distributions and, therefore, 
benefits from existing first-order techniques to simplify and expedite solutions with negligible consequences. 

B, Statistical Format 

In structural analysis, tolerance limit is a convenient statistical format for all input and output data 
which defines the distribution and specifies a statistical range of variations about the data’s most probable 
value. Statistical tolerance limits may be defined from a normal probability density plot for any given 
proportion of data and are commonly expressed in statistical parameter format: 

T Um =n±No , ( 6 ) 


or in using equation (4): 


T lim =H(l±Nn) , ( 7 ) 

where Tis the input-output data, such as loads or stress. The range factor N is a designer controlled number 
of standard deviations required to capture a specified percent of data from the high or low side (±) of a plot. 
On many nonredundant critical structures, the range factor of N— 3 is selected to capture a realistically high 
percent or observed extreme of the data population. In current deterministic analysis, the single values of 
the mean and standard deviation are substituted in equations (6) or (7) to reduce them to single value input- 
output data. It must be emphasized that once the tolerance limit format is converted into a deterministic 
value, it cannot be decomposed into its original statistical format from the deterministic value alone. 

However, true values of the mean and the standard deviation in equation (6) are not generally 
known from small sample sizes, for they may not contain a given portion of the population estimated by 
equations (2) and (3). In other words, the same test conducted on the same number of specimens by different 
experimenters will result in different means and standard deviations because of the inherent randomness in 
the specimens and testing. The population must contain results from all of these experiments. 

To insure with a certain percentage of confidence that the given portion is contained in the population, 
a A"- fact or is determined to account for the sample size and proportion. Figure 2 provides the A"- factor for 
random variables with 95 percent confidence levels and three probabilities (0.90, 0.95, and 0.99) in a one- 
sided normal distribution. Other confidence ^-factors may be computed from the program provided in 
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appendix B from which the A'-factor maximum or minimum design value may be determined for a specified 
probability and confidence. That allowable design value is the upper or lower tolerance limit defined by: 

Tnm=H±Ke . ( 8 ) 


5.5 

5.0 

4.5 
5 4.0 
re 3.5 

Li_ 

i 3.0 

2.5 

2.0 

1.5 


FIGURE 2. — /^-factors for one-sided normal distribution. 


C. Resistive Uniaxial Stress 

The resistive stress probabilistic distribution is a direct data characterization of material strength 
from a uniaxial stress test. While most of NASA material properties are specified by “A” and “B” bases, 
the tolerance limit of equation (8) is specified for the lower half of the distribution: 

Sr=Vr—K(7 % , (9) 

and by the number of test samples. Test samples may range from standard uniaxial tensile specimens 
through pressure bottles. It’s a trade between the initial cost of extensive material sample testing and the 
recurring cost of lost performance of global structures based on designing to a larger factor to compensate 
for small sample property predictions. 

Consistent with critical main structures and welds, the stress dispersion is often autonomously 
assumed as 3 a, or K= 3, requiring at least 32 test samples (figure 2) for an A-Basis material. An A-Basis 
property allows that 99 percent of materials produced will exceed the specified value with 95 percent 
confidence. The B-Basis allows 90 percent with the same 95 percent confidence. Figure 3 illustrates the 
probability and confidence plot for an A-Basis design. 



0 10 20 30 40 50 60 70 80 90 

Number of Samples, n 
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FIGURE 3. — One-sided normal distribution with A-Basis. 

Most normally distributed material properties are developed in tolerance limit format as in equation 
(9). However, they are more often reduced and published as deterministic single values, and cannot be 
decomposed again into tolerance limit format as required for reliability analyses. These published 
deterministic properties are a serious loss of existing but incompatible data for future applications of reliability 
methods. 


D. Combining Statistical Data 

Combining data that are statistically characterized variables and are mutually exclusive may be 
defined as a multivariable function by combining their dispersions through the following error propagation 
laws. 8 « 9 When two or more independent variables are added or subtracted, their means are added or subtracted 
and their standard deviations are root- sum -squared (rss) by the summation function rule: 


for z=x±y, (T z =^j<7^+<Ty . (10) 

Applying this to the sum of a string of tolerance limits gives: 


r=X7] = 2A*/+ 

(=i t=i 


/=! 


( 11 ) 


When independent variables are multiplied and/or divided (± exponent), their coefficients of variation are 
root-sum-squared according to the power function rule: 


for z=x a y b , ri z = Ja 2 ^+b 2 i]^ . (12) 

where rj- represents the coefficient of variation of that product. Elastic modulus and Poisson’s ratio are 
defined by multivariables having measured dispersions and must be combined by the power function rule. 
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Though these rules apply to all statistically formatted and manipulated input-output data through loads and 
stress computational processes, deviations are inaccessible for complying with these functional rules in 
deterministic methods. 


E. Interfering Distributions 

Another type of data characterization is the development of a third distribution from two opposing 
distributions defined by their tolerance limit formats, such as in the failure concept. Failure occurs when 
demand exceeds capability. When applied stresses and material strengths are defined by probability 
distributions, probability of failure increases as their tail overlap area increases as shown in figure 4. The 
overlap area suggests the probability that a weak material will encounter an excessively applied stress to 
cause failure. The probability of failure decreases as the designer controlled difference of the distribution 
means increases and the natural distribution shapes decrease. 



FIGURE 4. — Failure concept. 

As discussed before, just the distribution side producing the worst case design problem is of any 
engineering interest, as is clearly demonstrated by the failure concept of figure 4. Only data from the right 
half of the applied stress distribution (greatest demand) is engaged with data from only the left side (weakest 
capability) of the resistive stress. Data from the other two disengaged distribution halves is irrelevant to the 
failure concept. Having defined the resistive stress distribution in subsection 3 above, the applied stress 
distribution computational process follows before the failure, or reliability, concept may be developed in 
section V. 
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IV. APPLIED UNIAXIAL STRESS 


This section develops the normal distribution worst half of the uniaxial equivalent applied stress on 
a first-cut structural frame design, and then iterates allowable stresses to satisfy specified safety. It 
characterizes probability data in statistical format and scrutinizes all data input-output throughout structural 
computation processes for dispersions and assumptions contributing to excessive conservatism. Structural 
computational processes leading to applied uniaxial stress normal distribution include: multiaxial quasi- 
static design loads; multiaxial static design loads; multiaxial stresses; state of stress; failure criteria; applied 
uniaxial stresses. 


A. Vehicle Quasi-Static Design Loads 

Events that design vehicle substructures include on-pad assembly, liftoff, max Q, high-g, separation, 
etc. Launch vehicle forcing functions used to generate ascent generalized forces include: wind speed, 
shear, gust and direction; propulsion thrust rise, oscillations, and mismatch; thrust vector control angle and 
rate; vehicle acceleration and angle of attack; mass distribution; other special trajectory generated 
environments. 

Because input environments to response analysis are time-dependent and statistically characterized, 
the induced loads output is also time-dependent and of a statistical nature. The response histories at select 
grid points are illustrated in figure 5, in which a specific time event may produce a maximum internal load 
for a degree of freedom (DOF) at one grid point only. Other time events produce maximum loads at other 
grid points as shown. Where a maximum internal load response is identified at a grid point, the free-body 
diagram of the included substructure experiencing that maximum response is constructed with all time- 
consistent loads acting along the total system. This computational process for designing different parts 
through time-consistent and statistically dispersed loads is repeated for each substructure at each unique 
event time producing the maximum load response about each axis. The end product of the structural response 
to environmental excitations is a set of maximum design loads, or “limit loads,” and event times for all the 
system substructures and critical regions. 



FIGURE 5. — Time-dependent response. 


12 



Aerospace loads modeling uses established computational structural dynamics principles and solution 
techniques 10 ’ 1 1 for multi degrees-of-freedom (MDOF) structures. Models assume the structural system to 
be represented by a network of finite elements designated along the body possessing mass, damping, and 
stiffness. Natural and induced environments act as forcing functions at discrete grid points. The motion of 
the total structure is composed of a system of substructures which are expressed by the linear matrix 
differential equation: 


[M){X(t)}+[C]{x(t)}+[K]{X(t)}={F(t)} . 


The acceleration matrix, {X(7)}, is the time-dependent physical coordinates at DOF. The [M], [q, and [K\ 

coefficients are mass, damping, and stiffness matrices, respectively. Dispersion of input data to these 
coefficients is usually constant or negligible and may be defined as cr= 0 throughout the computations. The 
forcing function {F(t)} is a matrix of time-dependent environmental excitations acting along the structural 
body. Its columns represent time increments. Its matrix rows represent discrete grid points of body internal 
DOF at which natural or induced environments are acting at one instant of time. Data characteristics of 
these forcing functions are of a statistical nature, having notable dispersions and consequences. Subsequently, 
they must be defined in statistical format and must be further treated by error propagation laws of equations 
(10) and (12). 

However, these forcing functions are currently defined by deterministic variables and are bounded 
to be summed algebraically when combining through all the following major steps in the quasi-static 
computations. The combined summation includes the deviations which generate a major source of 
uncontrolled conservatism. 

The above matrix differential equation is comprised of a set of coupled equations of motion which 
may be uncoupled through the mode-superposition method to determine the response of a system to a set 
of forcing functions. The system’s undamped natural frequencies <»and mode shapes [0] are solved from 
the undamped eigenvalue problem, ([AT|-£t) 2 [M\) [<j>]= 0, to obtain the coordinate transformation: 

N 

1*^(0 , 

r-\ 


where q is the generalized coordinates for r= 1, 2, 3, ...N modes. The shape matrix [(f)] rows represent the 
mode shape values at each DOF grid point, and columns represent different mode shapes relating to each 
natural frequency. Substituting the coordinate transformation equation into the linear matrix differential 
equation and premultiplying by the transpose of the mode shape, results in the equation of motion in terms 
of modal matrices and generalized coordinates. Because of orthogonality, coefficient matrices are diagonal 
matrices, and the uncoupled system differential equation of motion reduces to: 


UM0}+[2 £co]{q(0}+[co 2 ]{q(t)}=[<t>] T {F(t)} , 
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where [/]=[0] 7 [M][0], [2£y H<pf [C][0], and [c o 2 ]=[<j)] T [K][<p ] are the generalized (unity) mass matrix, 
damping matrix, and stiffness matrix, respectively, in deterministic format, and [0] {F(t)} is the generalized 
force in statistical format. 

The generalized force is calculated from a given a set of forcing functions also in statistical format, 
and the generalized coordinates q,q,q are then determined by integrating the uncoupled system differential 
equation for events characterized by associated forcing functions. Substituting these generalized coordinates 
into the coordinate transformation equation above yields the desired system physical coordinates X,X,X. 
These physical coordinates are then used to compute the substructure internal loads to form a set of quasi- 
static design loads calculated through a loads transformation matrix (LTM) of the inquired internal loads. 
Applying the substructure’s stiffness matrix into the modal displacement method, the internal loads {L(t)} 
of the substructure are given by: 




where [/el selects rows of the substructure stiffness matrix corresponding to the desired internal DOF 
grid points, and: 

{X(t)} = [T]{X(t)} 

is the total substructure displacements. The [7] matrix selects the substructure DOF out of the system 
displacements. The loads matrix may be rewritten as: 

{L(t)} = [LTM]{q(t)} , (13) 

where: 

[LTM]=[fC I [71 [0] 

is the load transformation matrix. The resulting internal loads, {L(t)}, are the desired quasi-static response 
load at grid point “g” substructure, and results are expressed with forcing functions F, in statistical format 
and “c,” time consistent response gains (influence coefficients): 

L g =c \ F \ +C 2 F 2 +C 3 F 3 +C 4 F 4 +C 5 F 5 +K 
or: 

L =c, (/J { +N i CJ l )+C 2 (/J 2 +^2 (7 2 ^ +c 3^3 + ^3 (T 3^ +c 4^4 + ^4 (J 4^ +K 

o 

influencing respective subscript grid point. The resulting equation ( 1 5) defines a linear combination of the 
elements of a random vector having a combined mean: 


(14) 

(15) 


»8 = 




n 

i = 1 


(16) 
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and combined standard deviation: 


1 

2 




5W 

Z=1 


Kf'N'O'f 
i = I 


( 17 ) 


Combining equations (16) and (17), and autonomously selecting the final probability range factor N, the 
desired output probabilistic response (or limit) load at grid point g is: 

L g =n g +No g . (18) 

However, in the deterministic method, all the terms in equation (15) are summed which reduces the 
deterministic load response to: 


L g = lc i H i + lc i N i a i 
8 i=l /=! 1 1 1 


(19) 


While the first term on the right side of equation (19) may be reduced to the combined mean of equation 
(16), the second term violates the summation function rule, equation (11), and reflects the worse-on-worse 
input-output process which is excessively conservative. 

B. Total Combined Structural Loads 


Other structural loads that must be combined with the quasi-static loads are the static ground and 
flight environmentally induced loads which include thermal, pressure, vibration, acoustic, etc. These loads 
are statistically derived and must be statistically formatted as in equation (6) and combined consistently 
with operational event and time by equation (11): 


, ( 20 ) 

2 N k 

k=l 8 

or: 

L g =n g (l+N g rj g ) , (21) 

where L g is the load at grid point g, and k represents all the induced loads at an event and time. N g is 
autonomously selected for that grid point. Induced loads with negligible variations and consequences may 
simplify the computations by assuming <y= 0 in the summation function rule. 
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However, in the current deterministic method, all static loads are reduced to deterministic terms 
and added to deterministically reduced quasi-static loads in which all means and standard deviations are 
added: 

L r lU t +lN k <, k , (22) 

° £=1 k =\ 

contrary to the summation function rule of equation (11). Furthermore, the statistically formatted static 
load cannot be appropriately added to the deterministic quasi-static load. This step in the loads process 
generates a second source of excessive conservatism, in which the correct total load derived from equation 
(20) is smaller than the determinist total load from equation (22). 

C. Structural Stress Response 

The structural stress model is represented by the same network of finite elements designated along 
the body with grid points corresponding with those on the quasi-static loads model above. Then each 
substructure is analyzed for the worst combination of time constant loads acting over the grid points that 
produce the greatest applied stress. Since input loads are expressed in statistical format as in equation (21), 
the computed stress output should also be produced in statistical format: 

S=n g (l+N g ri g ) , (23) 

through the error propagation laws of equations (10) and (12) for all stress components at each grid point. 
The coefficient of variation t] g is the same as that derived from loads in equation (21). 

In deterministic methods, the input statistical variables, equation (21), are treated as single values, 
and they again violate the error propagation laws when combined in computations of the structural response. 
To create a third source of excessive conservatism the NESSUS/FEM 12 module should be considered 
when using finite element methods. 

D. Equivalent Uniaxial Strength 

Proceeding with the search for sources of excess conservatism in the current deterministic process, 
the state of stress and failure criteria were examined. In order that applied triaxial stresses acting at any grid 
point, may be equated to the resistive (or ultimate) uniaxial stress in equation (1), the applied triaxial 
stresses must first be reduced into one dimensional (resultant) stress and then indexed to an equivalent 
uniaxial yield strength. 

The complex state of stress at a point on an oblique surface of a solid may be readily derived^ by 
modeling the three normal principal stress components acting along the orthogonal principal axes of a 
tetrahedron. The sum of forces along each axis provides three linear homogeneous equations to be solved 
simultaneously. A nontrivial solution of stress on the oblique surface is obtained by setting the resulting 
determinant of the stress coefficients to zero. The solution to the determinant is reduced to a cubic equation 
having three combinations of component stresses as coefficients 7/ of the oblique normal stress: 
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s 3 =i ] s 2 -i 2 s-i 3 = 0 , 


known as invariant. The First invariant is the sum of the determinant diagonal which relates to the 
hydrostatic stress: 


7 j — 5 | +S2+S3 1 

with a mean stress of S mean =/|/3. The second invariant is the sum of the principal minors: 

I 2 =^[(S ] -S 2 ) 2 +(S 2 -S 3 ) 2 MS 3 -S x ) 1 } , 

that relates to shear stress. These shear stresses should be in statistical format but are currently defined as 
deterministic values. The third invariant is the determinant of the whole matrix. These invariants of the 
state of stress are noted to be independent of material properties, which incites the next point. 

Currently, there is no theory that directly relates multiaxial stresses with uniaxial yield or ultimate 
stress. However, there are several criteria in which the elastic limit of a multiaxial stress state is empirically 
related to the uniaxial tensile yielding, and results are reasonably consistent with experimental observations. 
The Mises yield criterion 14 is based on the minimum strain energy distortion theory which supposes that 
hydrostatic strain (change in volume) does not cause yielding, but changing shape (shear strain) does cause 
permanent deformation. Hence, the yield criterion relates the experimental uniaxial tensile elastic limit, S ty , 
to the principal shear stresses through the square root of only the second invariant of the stress matrix. Then 
using the second invariant, the Mises initiation of yield criterion is expressed in its familiar form by: 

• < 24 > 

which depends on a function of all three principal shear stresses. Because of squared terms, it is 
independent of stress signs and, therefore, it is applicable to compression and tensile combinations of 
multiaxial stresses. And because of isotropy, the second invariant implies that it is independent of 
selected axes and may be expressed about any oblique plane: 

^=[ s2+s2+s ?-V,-V ; -V2 +3 < s i + ^ +s2 ,T • (25) 
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Using equation (25), the pure shear yield stress reduces to S sy = -y=- and is a good approximation of test 
results. Having established the yield stress by equations (24) and (25), the criterion also expresses the 
equivalent applied uniaxial tensile stress over the total elastic and inelastic range about that yield stress: 


S ■ = 
eqi v. 




(26) 


As noted above, these invariant stresses are in fact probabilistic stresses. But the current deterministic 
application of the Mises criterion, as expressed by equation (26), again violates the error propagation laws 
enforcing a larger combined stress case than the statistically complied case, thus rendering a fourth source 
of excessive conservatism. Since the four sources of excessive conservatism build on each other sequentially 
in the applied stress computational process, equations (25) and (26) represent the total accumulative yield 
and allowable stresses respectively. If the Mises stress of equation (25) exceeds the allowable stress of 
equation (1), the structural thickness is increased and the applied stress process is iterated. 

Returning to the Mises criterion of equation (26), the local multiaxial stresses should be in statistical 

format: 


S) =(/!■+ N&) , 


(27) 


and may be appropriately combined through the error propagation laws by expanding the functional 
relationship in a multivariable Taylor series about a design point (mean) of the system. The mean of the 
Mises combined applied stresses is determined from equations (26) and (27): 


V A = tic +A ~H Z » X +3( \i% +fij z +/£ )] 2 • 


(28) 


The combined standard deviation is calculated from: 




\ 
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+ 
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s 2 
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as z 

V °°z . 
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as Gzx 

\ m zx J 


(29) 
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and the controlled standard deviation is: 
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V. SAFETY CHARACTERIZATION 


Having defined the probabilistic resistive stress by equation (9), and having amended the probabilistic 
applied stress by equation (35), a relative probabilistic safety factor may be established, and its experimental 
verification limits, its relative role in the failure concept and absolute safety index may be completely 
characterized. 


A. Probabilistic Safety Factor 

Generically, a safety factor is expressed as the ratio of the resistive to the applied stresses, where the 
resistive stress is reference to the ultimate or yield strength of an A- or B-based material and the applied 
stress is designed to not exceed an arbitrarily specified safety factor. A probabilistic safety factor expresses 
the applied and resistive stresses in probabilistic format: 


_V vJI-ktir) 

s a VaO+km 


(36) 


In designing to a specified safety factor, the designer controllable variables are the: material ^-factor by 
choosing the number of test specimens as discussed before; the applied stress range factor selected 
autonomously by the loads analysts; and maximum predicted applied stress to not exceed the specified 
safety factor. 

To verify the design safety factor of equation (36) on a substructure, two conditions must be 
experimentally confirmed, the maximum predicted operational applied stress and the structural response to 
the predicted applied stress. It must be noted that the maximum predicted operational applied stress can 
only be confirmed by field or flight tests, and it may be many flights before each unique event producing 
the maximum load environment at each different substructure is achieved. It should be further noted that if 
the design operational applied stress is exceeded well into the program operational phase flight (exceeds 
yield safety factor), the contingency safety factor of equation (1) will accommodate it as nonlinearly inelastic 
on the “first loading” (cycle) and subsequently linearly elastic. 

The structural response of the NASA standard is experimentally verified through a static test in 
which the maximum predicted operational applied stress with a safety factor of at least one will cause the 
material to yield, and the operational applied stress increased by a safety factor of at least 1.4 times will 
produce fracture. If yield and fracture occur prematurely, the stress math model or material properties are 
in error, or a sneak phenomenon may be involved. Response dispersions of test articles vary with materials 
and manufacturing dispersions. 

The probabilistic structural response of equation (36) is experimentally verified to a predicted leaner 
applied stress. Structures tested to the more conservative deterministic applied stresses possess more margins 
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than the more appropriate probabilistic applied stresses . Therefore, it is expected that those existing structures 
tested to deterministic stresses and found submarginal may be perfectly adequate. Those found adequate 
may be excessively conservative. 


B. Failure Concept 

While the prevailing deterministic safety factor is the most commonly specified structural safety 
criterion, the mechanics of how safety is achieved in the failure concept are not expansively understood. In 
applying the probabilistic safety factor to the failure concept diagram of figure 6, equation (36) is noted to 
define three zones: 


Vr-Sr = Vr-(Vr-Kv r )=K<7r , (37) 

S R -S A =(SF)S a -S a =(SF-\)S a =(SF-\)(fj. A +N a g a ) , (38) 

S A-^A=(^A +N A^A>-^A= N A< j A ■ (39) 


Applied 

Stress 

(Demand) 



Resistive 

Stress 

(Capability) 


FIGURE 6. — Failure concept governing zones. 


The mid zone is observed as an extension of the applied stress distribution by the safety factor: 

Sr — ^a = (^ ~ Sf)^R ' 

Thus, the primary role of the safety factor is to decrease the applied stress tail overlap in the failure concept 
by extending the applied stress tolerance limit in figure 6 through the combined equations (38) and (39): 

N eff° A = S R~Fa= N A°A +<SF-m A + N a O a ) , 

from which the effective applied stress range factor is: 



SF-1 

7 U 


+n a sf . 


(41) 
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The sum of the three zones is the difference of the applied and resistive stress distribution means: 

Hr-Va=(Vr-Sr)+( s r- S a) + ( S A-Va) • ( 42 ) 

Though equation (42) combines probability contributions of applied, resistive, and safety factor range 
factors which together decrease the tails overlap, and failure as shown in figure 6, their integrated difference 
suggests a relative safety assessment looking for an absolute safety reference. 

C. First Order Reliability 

Many advanced techniques are being investigated and are evolving 15 for providing reliable structural 
designs. In the meantime, assuming split normal probability distributions and defining resistive and applied 
stresses in probabilistic format as in equation (36) leads to the first order reliability method may be compatible 
with prevailing practices, codes, and skills. 

The concept of failure was introduced by figure 4, in which the distribution tails overlap suggests 
the probability that a weak resistive material will encounter an excessively applied stress to cause failure. 
This is to say that the probability of success is reliability and that the reliability is less than 100 percent. 
Therefore, the probability of interference is the probability of failure and is governed by the difference of 
their means, p R -p A - Increasing the difference of the means decreases the tail interference area. 

Given that the applied and resistive stress probability density functions are independent, they may 
be combined to form a third random variable density function 16 in y =Sr-Sa ■ If Sr and Sa are normally 
distributed random variables, then y =S-Sa are also normally distributed and: 



1 

cx,,V2/r 

y 


EXP 


1 

*1 

=c 

1 

1 

2 

<7 


L > JJ 


(43) 


where fl y = P R ~P A and <F y = yJc . The y-variable distribution is plotted in figure 7. 



FIGURE 7. — Density function of random variable y. 
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The reliability of the third density function expressed in terms of y is: 


R=P(S R >S A ) = P(y>0)=jP dy , (44) 

0 y 

y-fi v 

where Py is the y-density function of equation (43). Letting Z= — then o y dz=dy. The lower limit 
of Z is: y 


z ^l = Ejl Ha 
1 ° y 

As y approaches infinity, Z approaches infinity, and the reliability of equation (44) is reduced to: 


R=- 


v 


2k' 


Z 2 

ex P\ ~y 


\dZ 


(45) 


The integration of equation (45) is programmed in the appendix C. Given the reliability R, the 
safety index “Z” value is printed, which may then be translated into statistical design parameters through 
the safety index expression: 


z=-z,= . 


(46) 


Equation (46) formulates the probability concept. Increasing the safety index and the standard 
deviations increases the means difference, which decreases tail interference area and the probability of 
failure. The reliability relationship with the safety index is plotted in figure 8, and the reliability notation 
0.9„ represents n 9s after the decimal. 


0.9 2 

GC 0.94 

£ 

| 0.9 6 

2 

I 0.9 a 
0.9 10 

1 2 3 4 5 6 7 

Safety Index, 7 

FIGURE 8. — Reliability versus safety index. 
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Recognizing that the safety index from equation (46) shares the same difference and function of the 
applied and resistive stress distribution means expressed by equation (42), then substituting and simplifying 
reduces the safety index to the desired expression: 


(47) 


The designer autonomously selected variables are the safety factor and the range factors. The 
safety index is noted to be about an order of magnitude more sensitive 18 to the safety factor: 


than to the range factor: 




dZ 


^=C 6 N a t1a 


dN^ 
1 V, 


where the common coefficient is Q,= — — The safety index is even less sensitive to coefficients of 

Vr-Va 

variation. The coefficients of variation are passive variables derived from statistical data, and they vary and 
combine in the structural computational process from region to region with variations in applied stresses 
and materials (including welds), which suggests that the safety factor or safety index varies from region to 
region. A specified common safety factor will not provide a common reliability. 

The safety index provides structural analysts the options of designing substructures to specified 
uniform safety factors and then calculating the absolute reliability from region to region having different 
statistical variables, or of designing substructures to specified uniform reliability by adjusting safety factors 
from region to region. Equation (47) is particularly important to the structural system in that the combined 
probability of range and safety factors autonomously specified by the materials, loads, and stress disciplines 
may be integrated and optimized into a specified absolute reliability uniformly across the vehicle structure 
regardless of region to region material or operational environmental changes. 

Developing a reliability criterion is an on-going concern in the aerostructural community. Selecting 
an arbitrary standard reliability is no better than the current standard safety factor. It should be expected to 
be derived from some compelling physical or economic constraint, such as risk: 

risk = X C k , 

k 

where P* is the probability of the outcome and Q is the consequence of the outcome k . However, risk 
requires great efforts, skills, and data. The results of risk are too sensitive to simple changes in failure 
processes, assumptions, and perceptions. A component common reliability criteria development is more 
compatible with the general designer culture and current practices. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The cause and nature of the long-suspected excessive conservatism in the prevailing structural 
deterministic method have been identified as an inherent violation of the error propagation laws incurred 
when reducing statistical data to deterministic values and combining them through several structural 
computation processes. These errors are restricted to the applied loads and stress distribution computations, 
and because mean and variations of the tolerance limit format are added, the errors are positive, serially 
cumulative, and excessively conservative. Since the quasi-static structural deterministic conservatism varies 
with design and analysis and is generally indefinable, stress audits based on deterministic pass-fail safety 
factors are speculative, all of which should provide incentives for adapting probabilistic structural methods. 

While most probabilistic methods in development would circumvent errors of propagation laws, 
patching, and partially converting existing deterministic methods may be more expedient and would provide 
the familiarity, confidence, and correlation with current experience base desired in new approaches. 
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APPENDICES 


The following programs are presented for the structural analysts’ information, convenience, and 
library. Recognizing that programs are computer and software specific, the following are coded in Quick- 
Basic for their simplicity and application to pocket computers, and because of their easy conversion to 
other languages. This appendix is broken down as follows: 


Al. Normal probability density 
distribution program 

A2. Normality distribution test program 
A3. Normalizing skewed distribution 
(Split normal) 

B. A'-factor program 

C. Safety index programs 

D. Mises criterion program 


Al. Normal Probability Density Distribution 
Program 

4 NORMAL PROBABILITY DENSITY CURVE 

OPEN”CLIP:”FOR OUTPUT AS #1 

INPUT “MEAN =”,M 

INPUT “STD DEVIATION =”,SD 

INPUT “START =”,XS 

INPUT “FINISH =”,XF 

INPUT “INCREMENTS =”,NX 

DX=(XF-XS)/(NX- 1 ) 

FOR 1=1 TO NX 
X=XS + (I-I)*DX 
F=EXP(-.5*((X-M)/SD) A 2) 

F=F/((2*3.!4159*SD) A .5) 

WRITE #3,X,F 
PRINT X,F 
NEXT I 
CLOSE #1 
STOP 


A2. Normality Distribution Test Program 


4 Input data 
CLEAR: INPUT “N=”;N 
DIM A(N),D(N),Z(N) 
FOR 1=1 TON 
INPUT AG) 

NEXT I 


4 sort data 
K=N-1 

LINE180:FOR X=I TO K 
B=A(X) 

IF B<=A(X+1) GOTO Iine250 
A(X)=A(X+1) 

A(X+1)=B 

Y=1 

T=X-1 

Iine250:NEXT X 
IF Y=0 GOTO line300 
Y=0 
K=T 

GOTO LINE 180 
line300: 

PRINT “SORT DONE” 


4 NORMAL DISTRIBUTION TEST 
4 Kolmogorov-Smimov (normality test) 

4 Critical values (n > 30): a=. 10, d=.805; 

a=.05, d=.886; a=.01, d=1.03 
OPEN”CLIP:”FOR OUTPUT AS #2 


‘mean and std. deviation 
FOR 1=1 TON 
C=C+A(I) 
D=D+A(I)*A(I) 
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NEXT I 
M=C/N 

SD=((D-N*M*M)/(N- 1 )) A .5 
PRINT “MEAN=”;M 
PRINT “STD DEV=”;SD 

‘standardized normal 
FOR 1=1 TON 
Z(I)=(A(I)-M)/SD 
NEXT I 

‘cumulative normal 
FOR 1=1 TON 
X=Z(I):T=X 

G=EXP(-X*X/2)/SQR(2*3. 14159) 

A 1 =.3 1 938: A2=-.35656: A3= 1 .78 147 
A4=- 1.821 25: A5= 1.330427 
IF X<0 THEN T=-X 
Y=1/(1+.2316419*T) 

P=((((A5*Y+A4)*Y+A3)*Y+A2)*Y+A1)*Y 

F=1-G*P 

IF X<0 THEN F= 1 -F 

DI=I/N ‘empirical cumulative 

D(I)=ABS (F-DI) 

IF DM<D(I) THEN DM=D(I):J=I 
NEXT I 

‘results 

FOR 1=1 TON 
PRINT D(I) 

NEXT I 

PRINT “WORST SAMPLE #";J 
PRINT “ABS DIFFERENCE, D“;DM 
CLOSE #2 
END 

Figure A1 is a plot of the D critical values 
for a one-sided distribution. The distribution is not 
normal if the program test result exceeds the “D 
critical value. Most engineering data distributions 
are one-sided, occurring in the lower or upper sides. 



FIGURE A1 . — One-sided test critical value of D. 


A3. Normalizing Skewed Distribution (Split 
Normal) 

Stress data are assumed to be based on a 
series of observed measurements reduced into a 
frequency distribution, or probability histograms, 
shown in figure A2. The base of the histogram is 
bounded by successive and equal ranges of measured 
values, and the heights represent the number of 
observations (frequency) in each range. 



Stress Distribution, x. 

FIGURE A2. — Stress frequency distributions. 

To illustrate the direct normalization of a 
skewed distribution, the stress frequency distribution 
data of figure A2 is applied to equations (2) through 
(5). Because the greater stress side defines the worst 
demand case (applied stress), only data from the 
shaded right side is used to calculate the normalized 
distribution variables. The distribution may be 
normalized by constructing a mirror image of the 
engaged side about its peak frequency value and 
calculating the standard deviation from the 
constructed symmetrical distribution. 

The peak frequency from figure A2 distribution is 
the mean, jd =14 ksi. 

Sample size is £/!=-«] +2Xn ( 

= 8+2 (8+7+4+2+0=36. 
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Sum of variations about the mean 

v=2'Zn i (x i -/i) 2 
2x8(14.0-14.0) 2 =0 
2x7 (14.5-14.0)2=3.5 
2x4(15.0-14.0)2=8.0 
2x2(15.5-14.0)2=9.0 
2x1 (16.0-14.0)2=8.0 
v=28.5. 

The variance is 


<T 2 =^-r=%^ = 0.81 . 

2,« — l 35 

and the standard deviation from equation (3) 
is cr=0.90 ksi. 

The coefficient of variation from equation (40) 

is 77=^=^=0.065 . 

‘SPLIT NORMAL DISTRIBUTION 
‘INPUT DATA 

CLEAR: INPUT “PEAK FREQUENCY^ ’,MU 
INPUT “NUMBER OF BARS=’\N 
DIM F(N), X(N) 

FOR 1=1 TO N 
INPUT F(I) 

NEXT I 
FOR 1=1 TON 
INPUT X(I) 

NEXT I 

‘CACULATE 
S=0:FOR 1=2 TO N 
S=S+2 *F(I) 

NEXT I 
SS=S+F(I) 

PRINT’S AMPLE SIZE=’\SS 
VMU=0:FOR 1=1 TON 
VMU=VMU+2*F(I)*(X(I)-MU) A 2 
NEXT I 

PRINT’VARIATIONS FROM MEAN=”,VMU 

PRINT “MEAN=”,MU 

SD=(VMU/(SS- 1 )) A .5 

PRINT “STD DEV=’\SD 

COF=SD/MU 

PRINT’COEF OF VAR=”,COF 


B. /iT-Factor Program 

‘K- FACTOR 
MARIO: 

DEFDBL A-Z 

INPUT” SAMPLE SIZE=”;NS 

INPUT ‘ ‘PROPORTION” ;P 

INPUT “CONFIDENCE=”;CL 

IF NS >90 THEN PRINT’S AMPLE SIZE SHOULD 

BE SMALLER THAN 90":WHILE INKEY$=””:WEND 

START=TTMER 

PI=3. 141592654# 

‘INVERSE NORMAL 
Q= 1 -P:T=SQR(-2*LOG(Q)) 

A0=2.30753:A 1 =.2706 1 :B 1 =.99229:B2=.048 1 
NU=A0+A 1 *T:DE= 1 +B 1 *T+B2*T*T 
X=TNU/DE 

L0: Z=1/SQR(2*PI)*EXP(-X*X/2):IF X>2 GOTO L3 

V=25-13*X*X 

FOR N=I 1 TOO STEP- 1 

U=(2*N+ 1 )+(- 1 ) A (N+ 1 )*(N+ 1 )*X*X/V 

V=U:NEXT N 

F=.5-Z*X/V 

W=Q-F:GOTO L2 

L3:V=X+30 

FOR N=29 TO 1 STEP-1 

U=X+N/V 

V=U :NEXTN 

F=Z/V :W=Q-F :GOTO L2 

L2:L=L+1 

R=X:X=X-W/Z 

E=ABS(R-X) 

IF E>. 00001 GOTOLO 
‘END OF INVERSE NORMAL 

‘CALCULATION OF FACTORIAL 
N=NS:NU=N-1 

MT=INT(NU/2):UT=NU-2*MT 

GT=1 

FOR 1=1 TO MT- 1 +UT 
KT=I 

IF UT=0 GOTO LI 
KT=I-.5 
Li:GT=GT*KT 
NEXT I 

GT=GT*( 1 +UT*(SQR(PI)- 1 )) 

GF=GT*2 A (NU/2- 1 ) 

‘END OF FACTORIAL 

SECANT METHOD 
KP=X:J=1 :K=KP 

K0=K:GOSUB INTEGRATION :SF0=SF 
K=K*( 1 +.000 1 ):K 1 =K:GOSUB 
INTEGRATION : SF 1 =SF 


29 



BEGIN:K=K 1 -SF I *(K I -KO)/(SF 1 -SFO) 

IF ABS((K 1 -K)/K 1 )<.00000 1 GOTO RESULT 

J=J+1 :K0=K1:K1=K:SF0=SFI 

GOSUB INTEGRATION:SFI=SF:GOTO BEGIN 

RESULT:FINISH=TIMER 

BEEP:BEEP 

PRINT “K =”;USING”##.####”;K 

PRINT “TIME=”;FINISH-START;"SECONDS” 

‘END OF SECANT METHOD 

WHILE MOUSE(0)<>( 1 ):WEND 

GOTO MARIO 

INTEGRATIONS 1=0:L2=10 
IF N>40 THEN L2=20 
DL=KP*SQR(N):TP=K*SQR(N) 

Y=NU/2 

M=2:E=0:H=(L2-L I )/2 
X=L1 :GOSUB FUNCTION 
YO=Y:X=L2:GOSUB FUNCTION 
YN=Y:X=Ll+H:GOSUB FUNCTION 
U=Y:S=(Y0+YN+4*U)*H/3 
START:M=2*M 
D=S:H=H/2:E=E+U:U=0 
FOR 1=1 TOM/2 

X=L 1 +H*(2*I- 1 ):GOSUB FUNCTION 

U=U+Y 

NEXT I 

S=(Y0+YN+4*U+2*E)*H/3 

IF ABS((S-D)/D)>.00001# GOTO START 

SF=S/GF-CL 

RETURN 

‘END OF SIMPSON 

FUNCTION: Z=TP*X/SQR(NU)-DL 
T0=Z:G0= I /SQR(2*PI)*EXP(-Z*Z/2) 

A 1 =.3 1 938 1 5: A2=-.3565638: A3= 1 .78 1478: 

A4=- 1.821 256: A5= 1 .330274 
IF Z<0 THEN T0=-Z 
W=l/( I +.23 1 649*T0) 

PI =((((A5*W+A4)*W+A3)*W+A2)*W+A I )*W 

PH=1-G0*PI 

IF Z<0 THEN PH= 1 -PH 

Y=PH*X A (NU- 1 )*EXP(-X*X/2) 

RETURN 

C. Safety Index Programs 

‘SAFETY INDEX FROM RELIABILITY 

‘NORMIN (.5,P,1) 

DEFDBL A-Z 

LL: INPUT'Probabi I i ty =” ; P 

PI=3. 141593 


PI=3. 141593 

Q= 1 -P:T=SQR(-2*LOG(Q)) 

A0=2.30753:al =.27061 

B 1 =.99229:B2=.048 1 

NU=A0+al*T 

DE=1 +B 1 *T+B2*T*T 

X=T-NU/DE 

CUMULATIVE NORMAL 
LO: Z= I /SQR(2 *PI) * EXP(-X*X/2) 

IF X>2 GOTO LI 

V=25-13*X*X 

FOR N=ll TOO STEP- 1 

U=(2*N+ 1 )+(- 1 ) A (N+ 1 )*(N+ 1 )*X*X/V 

V=U:NEXT N 

F=.5-Z *X/V 

W=Q-F 

GOTO L2 

Ll:V=X+30 

FOR N=29 TO 1 STEP-1 
U=X+N/V 
V=U:NEXT N 
F=Z/V:W=Q-F:GOTO L2 

L2:L=L+1 

R=X:X=X-W/Z 

E=ABS(R-X) 

IF E> .001 GOTO LO 
PRINT “SAFETY INDEX IS” 

PRINT USING “##.####”;X 

GOTO LL 

END 

‘RELIABILITY FROM SAFETY INDEX 

‘NORMIN (0.5,P,1) 

DEFDBL A-Z 

‘TNPUT’P=”;P:PI=3 .141593 
PI=3.141593 

‘Q= 1 -P:T=SQR(-2*LOG(Q» 
‘A0=2.30753:A1 =.27061 
‘B 1 =.99229:B2=.048 1 
‘NU=A0+aI*T 
‘DE=I+B1*T+B2*T*T 
‘X=T-NU/DE 

‘CUMULATIVE NORMAL 
INPUT’X=”;X 

Z= 1 /SQR(2*PI)*EXP(-X*X/2) 

IF X>2 GOTO LI 

V=25-13*X*X 

FOR N=ll TOO STEP- 1 

U=(2*N+ 1 )+(- 1 ) A (N+ 1 )*(N+ 1 )*X*X/V 
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V=U:NEXT N 
F=.5-Z*X/V:F=1-F 
GOTO L2 
LI:V=X+30 

FOR N=29 TO 1 STEP- 1 
U=X+N/V 
V=U:NEXT N 
F=Z/V:F= 1 -F 
L2: 

PRINT F 
END 

D. Mises Criterion Program 

‘ERROR PROPAGATION METHOD; 

‘ MISES CRITERION 
DEFDBLA-Z 

INPUT’NUMBER OF NORMAL STRESS=’\NS 
DIM STATIC NSM(3),NSSD(3),NSNF(3), 
NSFD(3),LNS(3) 

FOR 1=1 TO NS 

PRINT “NORMAL LOAD MEAN(";I;”)=” 

INPUT NSM(I) 

PRINT'NORMAL LOAD STD. DEVIATION(“;I;”)=” 
INPUT NSSD(I) 

PRINT'NORMAL LOAD N-FACTOR(“;I;”)=” 
INPUT NSNF(I) 

NEXT I 

INPUT “NUMBER OF SHEAR STRESSES=”,MS 
DIM STATIC SSM(3),SSSD(3),SSNF(3),SSFD(3), 
LSS(3) 

FOR 1=1 TO MS 

PRINT “SHEAR LOAD MEAN(“;I;’>” 

INPUT SSM(I) 

PRINT “SHEAR LOAD STD. DEVIATION(“;I;’>” 
INPUT SSSD(I) 

PRINT “SHEAR LOAD N-FACTOR(“;I;”)=” 

INPUT SSNF(I) 

NEXT I 

‘CALCULATION OF SYSTEM MEAN 
S1=0:FOR 1=1 TO NS:S1=S1+NSM(I) A 2:NEXT I 
S2=0:FOR 1=1 TO MS:S2=S2+SSM(I) A 2:NEXT I 
MZ 1 = S 1 -NSM( 1 )*NSM(2)-NSM( 1 )*NSM(3) 
MZ2=MZ 1 -NSM(2)*NSM(3)+3*S2 
MZ= SQR(MZ2) 

‘CALCULATION OF DERIVATIVES 


NSFD( 1)=(2*NSM( I )-NSM(2)-NSM(3))/2/MZ 
NSFD(2)=(2*NSM(2)-NSM( 1 )-NSM(3))/2/MZ 
NSFD(3)=(2*NSM(3)-NSM( 1 )-NSM(2))/2/MZ 
FOR 1=1 TO MS:SSFD(I)=3*SSM(I)/MZ:NEXT I 
‘CALCULATION OF SUM OF SQUARES OF 
‘ NORMAL STRESSES 
S3=0:S4=0:FOR 1=1 TO NS 
S3=S3+(NSFD(I)*NSSD(I)) A 2 
S4=S4+(NSNF(I)*NSFD(I)*NSSD(I)) A 2 
NEXT I 

‘CALCULATION OF SUM OF SQUARES OF 
‘ SHEAR STRESSES 
S5=0:S6=0: FOR 1=1 TO MS 
S5=S5+(SSFD(I)*SSSD(I)) A 2 
S6=S6+(SSNF(I)*SSFD(I)*SSSD(I)) A 2 
NEXT I 

‘CALCULATION OF SYSTEM STANDARD 
‘AND EFFECTIVE DEVIATIONS 
SZ=SQR(S3+S5):SN=SQR(S4+S6) 

NE=SN/SZ 

‘CALCULATION OF SYSTEM COEFFICIENT 

‘ OF VARIATION 

ETA=SZ/MZ 

‘CALCULATION OF SYSTEM TOLERANCE LIMIT 
TL=MZ+(NE*SZ) 

‘CALCULATION OF MISES FUNCTION 
FOR 1=1 TO NS 

LNS(I)=(NSM(I)+NSNF(I)*NSSD(I)) A 2 
NEXT I 

FOR 1=1 TO MS 

LSS(I)=(SSM(I)+SSNF(I)*SSSD(I)) A 2 
NEXT I 

FM1=0:FOR 1=1 TONS 
FM 1 =FM 1 +LNS(I):NEXT I 
FM2=0:FOR 1=1 TO MS 
FM2=FM2+LSS(I):NEXT I 
FM= SQR(FM 1 +3 *FM2) 

PRINT “COMBINED APPLIED STRESSS =”;FM 

PRINT “MEAN =”;MZ 

PRINT “STANDARD DEVIATION =”SZ 

PRINT “EFFECTIVE N = ? ;NE 

PRINT “COEFFICIENT OF VARIATION = ’;ETA 

PRINT “TOLERANCE LIMIT =” ;TL 
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